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Abstract. The well-known Turing machine is an example of a theoret- 
ical digital computer, and it was the logical basis of constructing real 
electronic computers. In the present paper we propose an alternative, 
namely, by formalising arithmetic operations in the ordinary computing 
device, we attempt to go to the analytical procedure (for calculations). 
The method creates possibilities for solving nonlinear differential equa- 
tions and systems. Our theoretical computer model requires retaining 
a finite number of terms to represent numbers, and utilizes digit carry 
procedure. The solution is represented in the form of a segment of a 
series in the powers of the step size of the independent variable in the 
finite-difference scheme. The algorithm generates a schematic representa- 
tion that approximates the convergent finite-difference scheme, which, in 
turn, approximates the equation under consideration. The use of prob- 
abilistic methods allows us to average the recurrent calculations and 
exclude intermediate levels of computation. All the stages of formaliz- 
ing operations of the classical computer result in "the method of the 
computer analogy". The proposed method leads to an explicit analyti- 
cal representation of the solution. We present the general features of the 
algorithm which are illustrated by an example of solutions for a system 
of nonlinear equations. 
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1 Introduction 



Analytical and numerical methods complement each other. Although most of 
nonlinear equations have no analytical solutions, yet numerical methods allow 
us to obtain solutions using the computer. It is important to emphasize that 
the computer is a material object (really, the electronic computer). In other 
words, an idealized mathematical object needs, in most cases, a materialized 
mathematical scheme. The question arose: can we formalize the operations of this 
material object itself in such a way that the numerical (calculating) procedures 
will become analytical? 
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Due to the efforts of some outstanding mathematicians such as Alan Turing, 
the computer theory has a soHd base for constructing different types of comput- 
ing machines (see [1]). The recent general tendency is as follows: most investiga- 
tions aim at elaborating new computing devices, e.g. quantum computers [2] or 
developing analytical procedures (and languages) directed for improving com- 
puter software [3], i.e., their general purpose is to elaborate new logical computer 
schemes. Yet, we can discuss other proposals concerning analytical computing 
techniques (and we study not the logical but the arithmetical structure of the 
computer). 

We will consider a possibility to compute analytically using the basic prop- 
erties of the classical digital computer (our main results are presented in [4, 
5]). From our point of view, there are two general features of the computer 
when dealing with numbers: 1) numbers are represented as segments of a power 
series; 2) tlicac; is a digit carry procedure. These properties can be applied in 
the iterative or finite-difference scheme, where in the recurrence formulas the 
intermediate levels can be excluded. These properties allow us, for example, to 
represent solutions of differential equations (or systems of differential equations). 
At each level, we present the solution as the segment of a series in the powers of 
the step T of the independent variable (we focus on the finite-difference scheme 
approximating the differential equation under consideration). Thus one can con- 
struct an analytical procedure, i.e. "computing without a computer", if using 
the computer analogy makes possible to exclude the intermediate expression 
(as a computer, in fact, cxchidcs the intermediate values when operating with 
numbers). The construction of the proposed method does not only reduce the 
number of the arithmetical operations in calculations. It provides a solution in 
the explicit form (as a computer provides a solution in the numerical form, after 
executing many intermediate and "hidden" operations). The final analytic so- 
lution is treated as the limit when the small parameter r tends to zero. In this 
short paper we describe the general features of the proposed method. 



2 The T-computer model 

In the present paper we consider general characteristics of our approach which 

can be applied when calculations are needed. A new theoretical model of the 
digital computer can be proposed if there is a small parameter in the problem 
under consideration. In this case one can use the expansion of the unknown 
function in a series in powers of this parameter. The universal approach of the 
finite-difference schemes (for solving differential equations) possesses such a pa- 
rameter, namely a step on an independent variable. 

Consider the following Cauchy problem for the following ordinary differential 
equation: 

f = ./(y), yio)-yo, (i) 

where / (y) is an infinitely differentiable function with bounded derivatives. We 
will refer to the argument t as "time" . Consider a finite-difference explicit scheme 
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of the following form: 

Vn+i =yn + TG (t/„) , (2) 

where G is also infinitely difFcrentiable and is defined by the chosen mimerical 
method. We assume that t is small enough that scheme (2) is stable; hence, 
it converges to the solution of problem (1). One can see that after performing 
procedure (2) we obtain a segment of the power series (with powers of t), but 
the number of terms increases with n. We will search for a representation of the 
solution by this series with restricted number of terms. The digit carry procedure 
is introduced that guarantees that the absolute values of the coefficients of the 
segments of the series are less than 1/r. If the coefficient becomes greater or 
equal to 1/r, its part is carried to the more significant digit (i.e. left digit). If 
the absolute value of the coefficient is less than 1/r, we call it normalized; if not, 
then we call it not normalized and denote it with the tilde aj,„. Consider the 
following representation of t/„: 

p 

yn = ^Piai,nT\ (3) 
i=0 

where the factors /3j e {—1, 1} are chosen such that all Oj are not negative. This 

segment of the series provides a unique representation for the numbers with error 
O (r^). Suppose that we apply (2) to obtain j/„ from We cannot guarantee 

that for any i the inequality aj^„ < 1/r holds true. Thus, 

p 

i=0 

Next, we apply carry procedure starting from ap,„: 

ap,n = ap,n mod 

Spm = (apm - ap,n) T = flp^riT ^ ap,nT mod 1 = [ap,„r] . 
For other coefficients, we obtain the following: 

+ f^m07n+l^m + l.n) mod -, 

r 

The condition < < 1/r implies that in the segments of the monotony, all of 
the coefficients Oj increase monotonically; thus aj,„ < Ui^n+i- Consider function 
G in (2). Let us expand it into the Taylor series near y* = (ioao^n- 



oo ^ / p 

G {Vn) = G (y*) + E fc! I E l^rnam,nl 
1 ' \m=l 

Using (1.2) and (1.3), we obtain the following: 

Vn+l 



d'^Giy) 



dy^ 



y=y' 



P oo / p \ 

m=0 1 ■ \m=l / 



d'^G {y) 



dy'^ 



v=v 
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Thus a theoretical model of the proposed r-computer includes constructing an 
arbitrary positional numeral system, that is a mathematical auxiliary object that 
is useful in our approach; it deals with powers of the step of the argument. In 
contrast to the ordinary definition with the integer radix we consider the radix 
as the real number, namely, 1/r. Consequently, the r-computer is an imaginary 
device that represents the numbers in the r-numeral system and applies carry 
procedure to prevent overflow of the digits. 

3 The method of computer analogy 

The model of the theoretical computer will be completed if we will introduce 
important properties of the stochastic behavior of the less significant digits. 
It allows us to reduce intermediate operations. Namely, the sums of stochastic 
numbers will be changed by appropriate means (expected values) which will be 
theoretically evaluated. In other words, "the inverse Monte Carlo procedure" 
will be realized. Then we will have the full computer model and the approach 
called the "method of computer analogy" (CA). 

The numbers that are produced by a certain formula are often referred to as 
pseudo-random numbers. These numbers are widely used in a variety of simula- 
tions and in Monte Carlo methods in general. An ideal random number generator 
(RNG) should produce numbers with the expected value E[x] = 1 /2 and the 
variance Var (x) = 1/12. Consider one of the most popular modern algorithms 
for generating random numbers by a linear congruent RNG, see e.g. [6]: 

Xm = {hxm-1 + c) mod P, 

where &, c and P are integers and h,c < P. There are also more complicated 
formulae. Let us consider as a sample the following nonlinear system (and the 
initial problem): 

— =v^-u^ u{0) = l 

— =u^-2v, V (0) = 0. 

which has no an exact solution (in the ordinary sense). The finite-diS'erence 
method (2) if we use the first-order scheme gives the following: 

Un+l = + Tvl - Tul^, 1*0 = 1, 
Vn+l ^Vn+ Tul - 2tW„, Vq = 0. 

According to CA we search for a solution in the following form: 

It can be proved that for the quadratic nonlincarity it is sufficient to retain 
the third power in the segments of the series. It can also be proved that in the 
problem under consideration there is no carry to the zeros digits therefore the 
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first (linear term) digits have the main roles (see [4,5]). The carry procedure 
gives the following: 



Un+i = 1 - (ai,„ + 1 - 52,n) T + ^(a2,n + 2ai,„ - (53,„) mod J - 
- ^(as.n - &i,„ + ai,n + 2o2,n) mod , 

Vn+i = (6i,„ + 1 - W2,n) T - ({b2,n + 2ai,„ + 26i,„ - a;3,„) mod + 



(^{bs^n + ai,„ + 2a2,„ + 262,„) mod 



where 

^2,n = [(a2,n + 2ai,„ - 53,„) r] , uj2,n = [{b2,n + 2ai,„ + 26i,„ - a;3,„) r] , 

^3,n = [(^3,71 - + a?,„ + 2a2,„) r] , a;3,„ = [(63,„ + „ + 2a2,„ + 262,„) r] 

For simplicity denote a = ai and 6 = 6i. From this follows the expressions for 
the linear coefficients: 

n n 
m=0 m=0 

Here we suppose that ^2 and CO2 are random integer values with the probabilities 
governed by the coefficients a and b. 

This formula produces integer numbers in the segment [0, l/r — 1] . Fig. 1 
shows the behaviour of coefficients that depend on the number of the level (the 
discretized time). Coefficients ai,„ and &i.„ depend on ai,„_i and 6,:,„-i respec- 
tively, and they can also indirectly (due to carries) depend on the coefficients of 
the previous level. Coefficients ai^„ and 6i^„ do not exhibit stochastic behaviour, 
but they influence on all other values of the solution. Fig. 1 demonstrates the 
stochastic behaviour of coefficients aj,„ and The stochastic behaviour of 
02,71 and 62,71 slowly transform into regularity with oscillations. Nevertheless, we 
can still apply probabilistic methods when treating the carried values as weak 
dependent random values. 

The stochastic properties allow us to obtain the inverse dependence of the 
linear term on the discretized time. It appears that the coefficients of t in the 
nonlinear powers have a stochastic behaviour; thus, we can apply the statistical 
methods. The carried value ^2,n and L02,n are the only values that directly define 
the linear term. The next stage is expressing non-random terms rising the carried 
values. Instead of carried values, we use their expected values. We refer to these 
values as quasi-random numbers. Let us use the law of large numbers to average 
the sequence of these numbers using their expected values. The factor in the 
convergence formula for low discrepancy points may be contrasted with n~^/^ 
convergence (see [7]). So the analog of the law of large numbers gives the following 
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Fig. 1. Dependence of coefficients ai, 02, as, 61, 62, 63 on time. 



inequality: 

c 

< — , c = const, 
n 

Here ((5) denotes the expected value of a random variable 5. It is admissible 
to assume that convergence is ordinary rather than in probability. The numeric 
simulation confirms that the error is close to n~^. In Fig 2, one can see that the 
error does not exceed . plus some constant proportional to r. Thus the error 
can be estimated as c/n + O (r). Since r is of the order of we conclude that 
resulting error is of the order of n~^. 



E52,m E{52,m) 

m=l m=l 



4 Obtaining the explicit analytical form of the solution 

It is desirable to use the solutions without interpolation, therefore time intervals 
for u and v must be equal. For this purpose wc use a or 6 as a basis variable, 
then the other variable is scaled proportionally to the ratio of time intervals, 
because for a small time interval a and b are changed linearly. We will use a as 
a basis variable. 

Let us denote the expected values of 1 — ^2 and 1 — W2 as (a, 6) and 
{a, b). The time intervals where a and b are not changed are l/\Es (a, 6)| and 
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Fig. 2. Stochastic errors in CA. 



1/ \ Eu {a,b)\ respectively. The expected values are 



Eg (a, &) = (!- arf - b^r^, (a, b) = {1 - arY - 26t. 



From this follows that the expected values of 82 and 0^2 are constant on layers 
defined by a and b. Denoting the number of terms on a-th layer as Ana = 
rio+i — Ua, we obtain: 



a = g(l-^2,n)= 5^^(l-52,n). 



m=l 



Since the increment Aa equals 1, it follows that 



1 "i 2 

^5 (m) ^ (i) 
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As we have mentioned, we use a as an independent variable, thus n (the number 
of a level) and b must be expressed through a: 



a— 1 a— 1 



2t 



Ua = l-aT, i;„ = 6„T = T + T ^ J| i " (T ' I ' 



TO=1 k=m+l 

a-1 



(1 - krY 



1 

ta = Tna = T}_^ 



The explicit solution allows us to construct approximations for asymptotic cases. 
Such asymptotic approximation valid for small (and moderate) t is given in [5] . 
Namely 

a ^ 

Ua = l-aT, Va = baT = T + T{l- 2t {a + 1)) V , 

1 — 2mr 

m=2 

a— 1 ^ 
m=0 

The important issue is finding the limit of this solution when the step size tends 
to zero, it will provide the final form of the solution. The convergence of the 
solutions at different values of the step sizes is shown in Fig. 3 (in the top the 
scaled segment of v (t) is shown) . One can see that the accuracy of the method 
is of the order of t. 

The passage to the limit implies that the step size r tends to zero. Thus we 
obtain a limit form of the solution, which is a solution for the given problem in 
a certain range of the argument. Note, that the Sergeev computer [8], dealing 
with the infinitely small numbers, could be considered as a hypothetical tool for 
passing to the limit. 



5 Concluding renicirks 

In the present paper, we described the main features of a novel approach, that is 
based on the idea of a connection between analytical and computational math- 
ematics. A theoretical model of the computer was presented. We analyzed the 
possibility of formalizing the computer operations, and proposed a representa- 
tion of a solution in the form of a segment of the series in the powers of the step of 
the independent variable. This technique involves analytical work for obtaining 
the probabilities (which can be a difficult task). However, the main steps can be 
listed as follows: choose a convergent finite-difference scheme, find the number 
of terms to be retained, and obtain the expressions for the main coefficients in 
the representation of the solution. 

The proposed method of the computer analogy allows us to represent the 
solution of the problem as a convergent series that can be analyzed. 
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Fig. 3. Convergence to the solution. 
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